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Abstract. The famous Padberg-Rao separation algorithm for fe-matching 
polyhedra can be implemented to run in 0(|V| 2 |-E log(| V| 2 /|-E|)) time in the 
uncapacitated case, and in 0(| V\ \E\ 2 log(|V | 2 /|E|)) time in the capacitated 
case. We give a new and simple algorithm for the capacitated case which can 
be implemented to run in 0(| V\ 2 \E\ log(| V| 2 /|.E|)) time. 

Key Words: matching, polyhedra, separation. 



1. Introduction 

Let G — (V, E) be an undirected graph, let b £ be a vector of vertex capacities 
and let u € Z^J be a vector of edge capacities. A u-capacitated b-matching is a family 
of edges, possibly containing multiple copies, such that: 

• for each i € V, there are at most bi edges in the family incident on i; 

• at most u e copies of edge e are used. 

If we define for each edge e the integer variable x e , representing the number of times 
e appears in the matching, then the incidence vectors of u-capacitated 6-matchings 
are the solutions to: 

J2ees(t) x * ^ b i> for alii 6 V (1) 
< x e < u e , for all e e E (2) 
x e e Z, for all ee E. (3) 

Here, as usual, 8(i) represents the set of vertices incident on i. 

The convex hull in R E of solutions to (1) - (3) is called the u-capacitated b- 
matching polytope. Edmonds and Pulleyblank (see [Edm65] and [Pul73]) gave a 
complete linear description of this polytope. It is described by the degree inequalities 
(1), the bounds (2) and the following blossom inequalities: 



'eel. i ii i ^ r- ' 



x e + } Xf < 



b(W)+}Z feF u f 



2 

for all W C V, F C 6{W) with b(W) + ~£ feF u f odd. (4) 

Here, E(W) (respectively, S(W)) represents the set of edges with both end- vertices 
(respectively, exactly one end-vertex) in W, b{W) denotes y~2 ieW h. 

An important special case is where the upper bounds u e are not present (or, 
equivalently, Uij > max{bi,bj} for all {i,j} G E). The associated (uncapacitated) 
6-matching polytope is described by the degree inequalities, the non-negativity 
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Algorithm 1 Minimum T-cut [PR82] 
Input: 

Graph G, set T C V, and weights c G Qf . 
Output: 

A minimum T-cut. 

1: Compute a cut-tree for the graph G with weights c and terminal vertex set T. 

2: For each of the n — 1 edges of the cut-tree do 

3: Let 5(U) denote the cut induced by the cut-tree edge. 

4: Check the cut: 

Compute the parity |T (~l U\ mod 2 and the weight c(U) of the cut. 
5: If adequate, store U. 
6: End for 

7: Output the best T-cut U. 



inequalities x e > for all e G E, and the simplified blossom inequalities 



In their seminal paper, [PR82] devised a combinatorial, polynomial-time sepa- 
ration algorithm for 6-matching polytopes. A separation algorithm is a procedure 
which, given a rational vector x* G Q E lying outside of the polytope, finds a linear 
inequality which is valid for the polytope yet violated by x* . Clearly, testing if a 
degree inequality or bound is violated can be performed in linear time, so the main 
contribution of [PR82] is to identify violated blossom inequalities. 

For uncapacitated b-matching, Padberg & Rao reduce the separation problem to 
the computation of a minimum T-cut, for which they give a generic algorithm, see 
Algorithm 1. We will give the definition of the minimum T-cut problem in the next 
section. Abbreviating n := \V\ and m := \E\, this algorithm involves the solution 
of up to n — 1 maximum flow problems on a graph with n + 1 vertices and n + m 
edges. Using the well-known pre-flow push algorithm [GT88] to solve the max-flow 
problems, this leads to an overall running time of 0(n 2 mlog n2 /m). 

The Padberg-Rao separation algorithm for capacitated b-matching, however, is 
substantially more time-consuming. It involves the computation of a minimum T- 
cut on a special graph, the so-called split graph, which has up to n + m + 1 vertices 
and up to 2m + n edges. Up toji + m-1 maximum flow problems may be required 
to be computed. Using the pre-flow push algorithm, this leads to a worst-case 
running time of 0(m 3 logn). In 1987, [GH87] observed that the above-mentioned 
max-flow problems can in fact be carried out on graphs with only 0(n) vertices and 
0(m) edges. Although the idea behind this is simple, it reduces the overall running 
time for the capacitated case to 0(nm 2 log™ 2 /™). 

In this paper, we propose a new separation algorithm for the capacitated case 
whose running time is the same as that for the uncapacitated case. As well as being 
faster than the Padberg-Rao and Grotschel-Holland approaches, the new algorithm 
is much simpler and easier to implement. It also has a surprisingly simple proof of 
correctness. 

Our results also apply to the case of perfect capacitated b-matchings. 




for all W C V with b(W) odd. 



(5) 
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Algorithm 2 Blossom minimization 
Input: 

Graph G, set T C V, and weights c, d £ Qf . 
Output: 

A minimum blossom. 

1: Compute a cut-tree for G with weights min(c, d) and terminal vertex set V. 

2: For each of the n — 1 edges of the cut-tree do 

3: Let S(U) denote the cut induced by the cut-tree edge. 

4: Check the cut: 

Compute (3(U) as in (8). 
5: If adequate, store U along with the arg-min F. 
6: End for 

7: Output the best blossom (U,F). 



As well as being of interest in the context of matching, the algorithm has an 
important application to the Traveling Salesman Problem (TSP). The special blos- 
som inequalities obtained when hi — 2 for all i and u e = 1 for all e are valid for 
the TSP, and facet-inducing under mild conditions, see [GP79a], [GP79b]. Thus 
we obtain a faster exact separation algorithm for the TSP as a by-product. In fact, 
the algorithm is applicable to a general class of cutting planes for integer programs, 
called {0, i^j-Chvatal-Gomory cuts, see [CF96]. 

Parts of the contents of this paper appeared in the proceedings of the Xth IPCO 
conference [LRT04]. However, the proof of correctness of the algorithm is now 
substantially facilitated. 

2. Algorithms for minimum T-cut and blossom minimization 

Given a graph G — (V, E) , an even-cardinality set T C V and non- negative 
rational edge-capacities c £ Q+, the minimum T-cut problem asks for an odd cut 
(U, ClI) (where ClI is the complement of U in the vertex set) such that the set 
U C V is T-odd, i.e., \T D U\ is an odd number, and which minimizes, subject to 
this condition, the submodular function 

U ^ c{U) := c ?- 

eeS(U) 

In 1982, Padberg & Rao gave the first polynomial-time combinatorial algorithm for 
computing a minimum T-cut, see Algorithm 1. The key ingredient is the compu- 
tation of a Gomory-Hu cut-tree [GH61] in step 1. Given a graph G — (V,E), a set 
X C V, and non-negative rational vector of edge-capacities c £ Q^, a cut-tree with 
terminal vertex set X for G and c consists of a mapping n: V — > X with ir(x) = x 
for all x E X, and an adjacency relation ~ on the set X. (We adopt the conven- 
tion that the edges of G will be denoted by xy, and the edges of the cut-tree by 
x ~ y.) The adjacency relation shall make the set of terminal vertices into a tree. 
An additional condition is required to hold. Deleting an edge x ~ y of the cut-tree 
partitions the set X into two sets X x and X y , and thus defines a cut (U, U) in G 
by letting U := 7r _1 (A 2: ) and U := -K~ 1 (X y ). We call this the cut induced by the 
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edge x ~ y of the cut-tree. Now, the condition which is required is the following: 

for x, y € X with x ~ y, the cut induced by this edge of the cut-tree 
shall be a minimum (s, t)-cut in G with respect to the capacities c. 

With the algorithm given by Gomory & Hu, a cut-tree can be computed in time 
0{\X\nm\og^/ m ). 

In Algorithm 1, the time for "checking the cut" in step 4 is negligible (the values 
c(U) even come for free with the Gomory-Hu algorithm), and hence the Padberg- 
Rao method for computing a minimum T-cut runs in time O (\T\nm log n2 /m), as 
mentioned in the introduction. 

Now we come to the blossom separation algorithm of Padberg & Rao [PR82]. 
Reformulating and generalizing, we say that a blossom is a pair (£7, F) consisting 
of a set of vertices U C V and a set of edges F C 5(U) with the property that 
\T (~l U\ + \F\ is an odd number. Then, if two non-negative rational weight vectors 
c, c' G are given for the edges of G, the blossom separation problem is equivalent 
to the problem of producing a blossom whose value 



is strictly less than, one, if it exists. For the sake of completeness, we describe 
how this equivalence is established. Padberg & Rao [PR82] introduced, for each 
u G V, the term s u := b u — X^ee5(i) Xei w hich is the slack of the corresponding 
degree inequality computed with respect to a given vector x. Then they showed 
that the blossom inequality (4) can be re-written in the form: 



To decide if, for a given x, sets W and F exist which violate (7), we define, in a 
canonical and straight forward manner, a graph G* , capacities c and d and a set 
T of vertices of G* , in such a way that a blossom with value strictly less than one 
gives rise to a violated inequality (7) and vice-versa. Let G* be constructed by 
adding a new vertex v to G = (V, E) and connecting it with an edge vu to every 
u G V. Then for each e G E, we let 



For the edges vu of G*, we let c uv := s u and c' uv := oo. Finally, we define T as the 
set of all vertices u for which the value b u is odd, and we let v G T iff J2 U ^ s °dd- 
Now it is easy to see that for each blossom (U, F) in G* with v $ U, the inequality 
(7) with W := U n V is violated by 1 - /3(U, F). Note that 0(U, F) = /3(CU, F). 

As mentioned above, the blossom separation Algorithm of Padberg & Rao [PR82] 
is considerably more complex than the minimum T-cut algorithm. It requires to 
construct a special graph G with m + n vertices and 2m edges, on which then a 
minimum T-cut is computed. 

We now give an algorithm for what we call the blossom minimization problem: 
given G, T and c, c' as above, find a blossom (U, F) which minimizes (3(U, F). The 
blossom minimization algorithm is displayed as Algorithm 2. 



P(U,F):= ]T c e + E4 



eeS(U)\F feF 
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For fixed U C V, it has been observed by Padberg & Rinaldi [PR90] that 

/3(f7) :=min{/?(£/,F) F c S(U), \T n U\ + \F\ odd} (8) 

can be computed in time 0(\S(U)\) by first tentatively taking F := {e G S(U) \ 
c' e < c e }. Now if T n U\ + \F\ is odd, we have found a minimizing F. Otherwise, 
find / G 6(U) minimizing |c/ — c'j\ over / G S(U), because then the symmetric 
difference of F and {/} minimizes (3(U,-). 

This implies that the loop 2-6 in Algorithm 2 runs in time 0(n 2 ) and that the 
running time of Algorithm 2 is dominated by the computation of the cut-tree in 
step 1, which amounts to 0(n 2 m log" 2 /" 1 )- 

The similarity between the Padberg-Rao minimum T-cut Algorithm 1 and our 
blossom minimization Algorithm 2 is striking. Moreover, in the next section, we 
give a short and elegant proof of correctness of Algorithm 2, which is similar to a 
proof of correctness of Algorithm 1 given by Rizzi [Riz02]. At this point, we might 
note that /?(■), unlike c(-), is not in general submodular. 

3. A SIMPLE PROOF OF THE CORRECTNESS OF ALGORITHM 2 

Let a cut-tree for G with terminal vertex set X be given, where X D T. We 
say that an edge x ~ y of the cut-tree is T-odd, if the sets of the bipartition of 
X defined by x ~ y are T-odd. Thus, the set of T-odd edges of the cut-tree form 
what is called a T-join, and an edge in the cut-tree induces a T-cut in G if and 
only if the edge is T-odd. The next theorem is the keystone of the correctness of 
Algorithm 1. For the sake of clarity, we repeat the proof of [Riz02]. 

Theorem 3.1 ([PR82]). One of the T-odd edges of the cut-tree induces a minimum 
T-cut in G. 

Proof. Let U be a minimum T-cut. Now U is a T-odd set, hence there exists an 
odd number of T-odd cut-tree edges leaving T n U . Let x ~ y be one of them, and 
let S be the minimum (x, y)-cut it induces by (6). Since U is an (x, y)-cut, we have 
c(S) < c(U), and since x ~ y is an T-odd edge, S defines a minimum T-cut. □ 

Now we come to the proof of correctness of Algorithm 2. 

Theorem 3.2. One of the edges of of the cut-tree computed in Algorithm 2 induces 
the a set U which minimizes /?(•). 

Proof. Let U be a set which minimizes (3(-). Further, define the set T' as the 
symmetric difference of T with all sets {u, v} for all e = uv G E c' e < c e . 

Case 1: U is T'-odd. The proof of Theorem 3.1 shows that there exists a T'-odd 
edge of the cut-tree which induces a minimizcr of /?(•). 

Case 2: U is not T'-odd. Let / = x'y' G S(U) have the minimal value of |c/ — c'j\ 
among all edges in S(U). On the path from x' to y' in the cut-tree, at least one edge 
x <~ y has one end in U and the other not in U. Let S be the minimum (x, y)-cut 
defined by this edge. Abbreviating w := min(c, c'), we then have 

0(U) = w{U) + \c f - c' f \ > w(S) + \c f - 4| > f3(S). 

The first inequality holds since U is an (x, y)-cut. As for the second, if S is T'-odd, 
then S minimizes [3 since w(S) < w(S) + \cf — c' f \ < (3(U); but if |T' n S\ is even, 
then (S, {/}) is a blossom whence w(S) + \c f -d f \ = 0{S, {/}) > 0(S). □ 
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